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g ABSTRACT 

o 

We present an explanation for the unusual temporal feature of the GRB 030226 
afterglow. The R-band afterglow of this burst faded as ~ t 2 in ~ 0.2 days after 
the burst, rebrightened during the period of ~ 0.2 — 0.5 days, and then declined with 
, ~ t~ 20 . To fit such a light curve, we consider an ultrarelativistic jetted blast wave 

expanding in a density-jump medium. The interaction of the blast wave with a large 
density jump produces relativistic reverse and forward shocks. In this model, the 
observed rebrightening is due to emissions from these newly forming shocks, and the 
^ ' late-time afterglow is caused by sideways expansion of the jet. Our fitting implies that 

the progenitor star of GRB 030226 could have produced a stellar wind with a large 



density jump prior to the GRB onset. 
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■ 1. Introduction 

The gamma-ray burst (GRB) 030226 was detected by the High Energy Transient Explorer 
2 Satellite {HETE-2) at 03:46:31.99 UT on 26 February 2003 (Suzuki et al. 2003). The burst 
duration in the 30-400 keV band was beyond 100 seconds, and the fluence in this energy band 
was ~ 5.7 x 10~ 6 ergs cm -2 . Owing to the rapid localization of HETE-2, an early afterglow of 
this burst was observed. Fox, Chen &: Price (2003) detected an optical counterpart ~ 0.11 days 
after the burst, and Price, Fox & Chen (2003) found about 9 minutes later that it had faded. The 
subsequent spectroscopic observations revealed some absorption lines, indicating that this burst 
was at a redshift of z = 1.986 ± 0.001 (Greiner et al. 2003a). Thus, the isotropic 7-ray energy 
release is estimated to be ~ 6 x 10 52 ergs, assuming a cosmological model with Q m = 0.3, = 0.7 
and Ho = 65kms _1 Mpc -1 . 

The afterglow light curve of GRB 030226 has one unusual feature. The R-band afterglow 
declined as ~ t~ 2 shown by the four observations during the period of ~ 0.13 — 0.20 days after 
the burst, and about seven hours later, it began to decay with another simple power-law of t -2 
as late-time afterglows from jets usually do. However, contrary to jetted afterglows observed 
previously, the extrapolation of the early afterglow light curve of GRB 030226 crosses with the 
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late-time light curve at ~ 2 days rather than ~ 0.5 days, which implies an obvious rebrightening. 
Such a feature could result from post-burst energy injection. A few variants of this model have 
been studied (e.g., Dai & Lu 1998a; Panaitescu, Meszaros & Rees 1998; Zhang & Meszaros 2002; 
Nakar, Piran & Granot 2003). In this Letter, we explain the GRB 030226 afterglow light curve 
by assuming this burst to be in a density-jump medium (see §2). In §3 we discuss scenarios for 
density-jump formation and give order-of-magnitude estimates of the density-jump factor, which 
are consistent with required by our fitting. Our results are summarized in §4. 

2. The Model 

Before proposing our model, we first discuss implications of the observations on the GRB 
030226 afterglow: The observed light curve index (a^) after ~ 0.5 days since the HETE-2 trigger 
suggests that this afterglow might have originated from a highly collimated relativistic shock, 
because the transition of this shock from the spherical to spreading phase could lead to steepening 
of the light curve to the flux proportional to t~' p for p > 2, where p is the spectral index of the 
shock-accelerated electrons (Rhoads 1999; Sari, Piran & Halpern 1999; for the case of 1 < p < 2 
see Dai & Cheng 2001). Thus, p = 02 — 2 for GRB 030226. If the early afterglow is assumed to 
result from the shock expanding in a medium with density of n w oc R~ k , then the light curve index 
in the slow-cooling regime becomes ol\ = 3(p — l)/4 + fe/(8 — 2k) (Dai & Lu 1998b; Meszaros, Rees 
& Wijers 1998). The observed value of ol\ ~ 1.2, combined with the inferred value of p, implies 
k ~ 2, revealing a wind surrounding the burst. This environment as an expected signature of 
massive progenitor stars of GRBs was discussed theoretically in detail (Dai & Lu 1998b; Meszaros 
et al. 1998; Chevalier Sz Li 1999, 2000) and confirmed observationally for some bursts (see, e.g., 
Price et al. 2002 for GRB 011121). Recently, the wind interaction model was also shown to 
explain well the early afterglow of GRB 021004 (Li & Chevalier 2003). 

Therefore, we consider an ultrarelativistic jet with an isotropic energy of E, which first 
expands in a circumburst wind of the density profile of n w (R) = AR~ 2 with A = 3 x 10 35 A* cm -1 
for R < Rq. Here A* = (M/10~ 5 M Q yr _1 )/(T4,/10 3 kms _1 ), where M and V w are the mass loss 
rate and velocity of the wind respectively. After the internal shock emission, the jet will start to 
sweep up its ambient wind, leading to an ultrarelativistic blast wave. The initial hydrodynamics 
and emission were recently discussed by Wu et al. (2003) and Kobayashi & Zhang (2003). We 
further assume this blast wave which will hit an outer high-density homogeneous region for 
R > Rq. Dai & Lu (2002, hereafter DL02) have analyzed the hydrodynamics and afterglow 
emission from a relativistic blast wave in detail when it expands in a density-jump medium, which 
might have been produced around ~ 10 18 cm due to the interaction of a stellar wind with the 
ambient matter (see §3). 

The interaction of a relativistic blast wave with a large density jump can be described through 
two shocks: a reverse shock that propagates into the hot shell (viz., the R < Rq wind matter 
swept up by the blast wave), and a forward shock that propagates into the outer high-density 
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medium. Thus, there are four regions separated in the system of interest by the two shocks: 
(1) the unshocked high-density medium, (2) the forward-shocked high-density medium, (3) the 
reverse-shocked hot shell, and (4) the unshocked hot shell. We denote 7j as the Lorentz factor of 
region "i" measured in the local medium's rest frame, and n\ is the baryon number density of the 
medium at R > Rq (region 1). Defining the ratio (/) of the energy densities of regions 4 and 1, 
DL02 have got 

72 = 73 = 3l/4 , (1) 

and found that a density jump with a factor of much more than 21 can lead to a relativistic 
reverse shock. In this case, the Lorentz factor of regions 2 and 3 can further be scaled as 
72 = 73 = 1.4(1 + zy^Al^n^S^^iR/Ro)- 1 , where m, 3 = rai/10 3 cm~ 3 , and t is the 
observed time (in 1 day) at R = R = 8.0 x 10 17 (1 + z^^E^ 2 A~ 1/2 tl /2 cm (where E 53 is 
the isotropic-equivalent energy of the jet in units of 10 53 ergs). Comparing 72 (or 73) with 
the Lorentz factor of region 4 at the crossing time, 74 = 8.8(1 + z) 1 / 4 E 1 5 ^A~ 1/ % 1/4 , it is 
easy to see that the Lorentz factor of region 3 is much less than that of region 4 for typical 
parameters, showing that most of the initial kinetic energy is converted into thermal energy 
by the shocks. The radius at which the reverse shock has just crossed region 4 is derived as 
R A = R [l + 1.25 x 10~ 2 (1 + zf^E^ 12 A^l 12 ^ 1 ' 2 ] (DL02). Thus, R A /Ro - K 1 for typical 
parameters, so that the Lorentz factor of region 2 at R = i?A is well approximated by its Lorentz 
factor at R = Rq. If the initial half opening angle of the jet is 9q ~ 0.1 as in most of the afterglows 
analyzed by Frail et al. (2001), then it is possible that after the crossing time the Lorentz factor 
of the jet is far below Oq 1 and therefore the flux decays as t~ p due to sideways expansion as long 
as the jet is still relativistic (Rhoads 1999; Sari et al. 1999). 

In the standard afterglow model (Sari, Piran & Narayan 1998), the synchrotron flux 
F U; i of region "i" in optical to X-ray bands is determined by three quantities: the typical 
synchrotron frequency f m ,i) the cooling frequency v c ^ and the peak flux F VjTnaXt i . We divide 
the whole evolution into three stages: (I) before, (II) during, and (III) after the crossing 
time. Using Chevalier & Li's (2000) results, we get these three quantities at stage I: v^it) = 

4.8 X 10 13 (1 + z) 1 /2 e 2_ ie V 2 _ i ^^l/2 r 3/2 Hz5 u l {t) = g g x 1Q 13 (1 + ^-3/2^3/2^1/2^-2^/2 Rz 

and Fl jmax (t) = 0.59(1 + zf^e^E^A^Dl 2 ^ Jy, where £ = (p - 2)/(p - 1), e e = 0.1e e ,_i 
and 6b = O.le^-i are constant fractions of the internal energy density going into the electrons 
and the magnetic fields respectively, t is the observer's time in units of 1 day, and Dl^s is the 
luminosity distance to the source in units of 10 28 cm. 

At stage II, the spectral break frequencies and the peak flux evolve as 

1.3 x 10 12 (1 + z)el_ l e 1 I i 2 _ 1 C 2 A^ /2 tQ 2 {R/R )- 4 R Z , if i = 2, 
<, = <! 1.0 xl0 15 e 2 _ 1 e^' i e 2 i?53^ 1 n^ 1 (i?/ J Ro)- 2 Hz, if i = 3, (2) 
u l m (t )(R/R y A , if* = 4. 
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F 11 ={ 

u,max,i 



_ f 1.0 x 10 12 (1 + z)- 1 e B 3 Z 2 1 ^ 1 n 1 j /2 (t/t )~ 2 (i?/^o) 4 Hz, if i = 2,3, 
"\ ^(t )(i?/^o)- 4 , if » = 4. 

470(1 + z)V24/ 2 _ i ^V2^ 1 n l! 3^ /2 J D£ 2 28 [(i?/ J Ro) 3 - l](^o)- 2 Jy, if * = 2, 
28(1 + z)e 1 ^ 1 E 53 n 1 1 /2 DJ j 2 28 [(R/Ro) 2 - l](i2/12o)- 2 Jy, if t = 3, (4) 

. Fl, max (t K(R/R y 2 , if i = 4, 

where ( is the ratio of the total electron numbers of region 4 at radius i? G (i?o, -Ra) and at radius 
Rq. It should be emphasized that 1^4 in equation (3) is in fact a cutoff frequency because of 
adiabatic expansion of region 4, implying that there is no emission from this region at frequencies 
above 1/* 4 

From DL02, we obtain the time at which R = R& as follows: iA,2 = *A,3 = 2.0io for regions 2 
and 3. If the afterglow originates from a jet, one expects that at stage III the break frequencies 
and the peak flux of region 2 are given by v™ 2 = ^^(^a^K^Aa^) -2 , ^2 = ^2(^,2), and 
^max,2 = -^,max,2(*A,2)(*/iA,2) _1 (Sari et al. 1999). For the parameters adopted below, the 
cutoff frequency of region 3 is less than the optical frequency so that there is no emission from this 
region at frequencies above the optical band. 

Figure 1 presents an example case of our fitting to the multiwavelength light curves of the 
GRB 030226 afterglow. Please note that an error of 0.05 mag for the initial four afterglow data 
points in the figure has been taken as the error in the observed brightness of two calibration stars 
by Garnavich, von Braun & Stanek (2003). This error is smaller than a conservative estimate of 
0.3 mag by Zeh et al. (2003). The model parameters are taken: e e -\ = 0.9, €b-i = 4, £53 = 1.0, 
A* = 0.01, p = 2.01, to = 0.22 days, and = 0.1. If the emission off the line of sight (LOS) 
were neglected, one would see an initial drop at t ~ to, a rapid brightening at to <^ t <J £a,2 and an 
abrupt decay at t ~ £a,2 for all the light curves (DL02). However, such features are expected to be 
smoothed by the off-LOS emission. In Figure 1, we have considered the off-LOS emission effect for 
region 4 due to its ultrarelativistic expansion (e.g., 74 ~ 53 3> 1 for the taken model parameters) 
at t ~ to, following Kumar & Panaitescu (2000). However, this effect is insignificant for region 3 
because 73 ~ 2.9 is not far greater than unity at t ~ iA,2- 

We next discuss implications of the fitting. First, because the late-time afterglow at 
high frequencies in our model results from fast-cooling electrons, the spectral index should be 
(3 = p/2 ~ 1.0, which is consistent with the observed values, /3 pt,ob = 0.99 ± 0.42 at the optical 
band (Nysewander et al. 2003) and (3x,ob = 1-0 ± 0.1 at the X-ray band (Pedersen et al. 2003; 
Sako &: Fox 2003), respectively. Second, the jet's isotropic-equivalent energy (~ 10 53 ergs) from 
our fitting is close to the isotropic 7-ray energy release (~ 6 x 10 52 ergs). This shows that the 
efficiency of 7-ray emission from internal shocks is as high as ~ 40%. Third, our fitting gives a 
weak wind A* ~ 0.01, which implies non-detection of an optical flash from GRB 030226, as in most 
of the GRBs observed previously (Wu et al. 2003). Finally, the required value of ri\ indicates that 
the progenitor star of GRB 030226 could have been in a dense medium such as a giant molecular 
cloud (GMC). 
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3. Density-Jump Formation 

In this section we show that a large density jump seems to form in the vicinity of a massive 
progenitor star. A natural scenario for density-jump formation is the interaction of a wind from 
a massive star with its outer environment. This interaction produces a stellar wind bubble. 
Numerical simulations of the bubble dynamics in a homogeneous interstellar medium were 
performed by Ramirez- Ruiz et al. (2001), and here we give analytical studies to obtain the 
order-of-magnitude factor of a density jump by considering two kinds of environment. First, we 
discuss a simple case in which the outer environment of a massive star is a homogeneous GMC 
with density n\. Owing to the presence of two shocks (viz., a reverse shock that propagates into 
the wind gas and a forward shock that propagates into the GMC gas), the system has a four-zone 
structure consisting of (1) the unshocked GMC gas, (2) the shocked GMC shell, (3) the shocked 
wind gas, and (4) the unshocked wind gas moving with a velocity V w and with a mass loss rate M, 
where zones 2 and 3 are separated by a contact discontinuity. Assuming an adiabatic expansion 
and two strong shocks (Castor, McCray & Weaver 1975; Wijers 2001; Mirabal et al. 2003), we get 
the forward shock radius 

i4 ~ 4.8 x 10 18 Ml/V^Vj /5 <5 cm, (5) 

and the reverse shock radius 

^.r.OxlO^V^f^^ (6) 

where M_ 7 = M/1O~ 7 M yr -1 , V w ^ = T4,/10 3 kms _1 , m j2 = ni/10 2 cm~ 3 , and t w>5 is the 
wind lifetime in units of 10 5 years. Here the reason for scaling the stellar mass loss rate with 
1O~ 7 M yr 1 is that A* ~ 0.01 in our fitting implies M ~ 10 7 V w ^Mq yr 1 . This scaling is typical 
for main-sequence stellar winds, and further not far below typical mass loss rates of Wolf-Rayet 
stars (Willis 1991). One large density jump appears at the contact discontinuity, whose factor is 
approximated by 

X = -^TT - 1-7 x l^MZrV^ttll- ( 7 ) 

The radius at which the density jump is required to appear in our fitting is estimated as 
Rq ~ 2.2 x 10 18 cm, which is larger than R r sh and closer to R f sh for reasonable wind parameters, 
so the afterglow shock wave would be able to reach the density jump at the early time when the 
rebrightening was observed in the GRB 030226 afterglow light curve. 

Second, the surrounding medium of a fast wind from a massive star (e.g., a Wolf-Rayet 
star) may be a dense slow wind, because this star could have passed through a red supergiant 
(RSG) phase. During such a phase, the star has produced a wind with a larger mass loss 
rate Mrsg and a slower velocity Vase- This wind has formed a medium with density 
riRsciR) = 3 x lO^MRsc.W^RSG.i^cm- 1 , where Mrsg,-5 = M RSG /1O- 5 M yr" 1 
and Vrsg.i = ^RSc/lOkms -1 . The interaction of a RSG wind with its homogeneous 
environment (e.g., a GMC) could have given rise to a RSG wind bubble. The dynamics of this 
bubble is similar to the case discussed above, provided that they are both adiabatic. From 
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equation (5), we can write directly the radius of a thin shell swept up by the RSG wind as 
-Rsh.RSG — 1-9 x 10 18 M^gQ _ 5 Vrsq i^i^^^rsg 5 cm ' wnere *RSG,5 is the RSG wind lifetime in units 
of 10 5 years. Also, we obtain the factor of a density jump appearing at the contact discontinuity 
in the RSG wind bubble, xrsg ~ ii^RSG^-s^RSG.i^i^RSG.s- 

Furthermore, a fast wind from the Wolf-Rayet phase sweeps up the pre-existing slow wind, 
leading possibly to another wind bubble. Chevalier & Imamura (1983) gave self-similar solutions 
of this bubble. Assuming an adiabatic evolution and two strong shocks (Garcia-Segura & Mac 
Low 1995), we have the forward shock radius 

R { sh ~ 8.6 x W 18 M% 3 V^M-^V^% A t w , 5 cm, (8) 

and the reverse shock radius 

R r sh ~ 1.6 x 10 1& M^V^M^^^ cm - ( 9 ) 

Letting Rl h in equation (8) be equal to .R s h,RSG> we define a critical lifetime of the Wolf-Rayet 
wind, 

C - 2.2 x l^MZi^ 3 ]^ t _ 5 V^ tl nT^^ yrs. (10) 

If tw ^ twi the Wolf-Rayet wind only sweeps up the RSG wind. In this subcase, a large 
density jump occurs at the contact discontinuity and its factor is \ w — 7ZRgG(i? s j 1 )/7i J1 )(-Rgj 1 ) ~ 
360M_7^ 3 V^ 3 3 M^gQ _5^rsgi' wnere -^sh an d -^sh are given by equations (8) and (9) respectively. 
When R[ h ~ R sh . 

RSG; the above two density jumps merge to one larger density jump with a factor 

of 

X ~ XwXRSG ~ 4 x 10 3 MI 7 2/ V^^ ) _ 5 ^ 5 il n?54 5 Gj5 . (11) 

On the other hand, if t w > t™, the Wolf-Rayet wind not only sweeps up the RSG wind but also 
the outer GMC gas, in which subcase the dynamics can be approximated by equations (5) and (6), 
and thus the factor of a large density jump that occurs at the contact discontinuity is estimated 
by equation (7). Because Rq is larger than R r sh and closer to i?^ h in both subcases, the afterglow 
shock wave would be able to arrive at the large density jump at an expected early time of GRB 
030226. 

4. Conclusions 

In this paper we have shown that density jumps (or bumps) are likely to occur in the vicinity 
of a GRB, especially if a stellar wind has interacted with the ambient matter or if winds ejected 
at different velocities have collided each other, prior to the GRB onset. We have considered an 
ultrarelativistic jetted blast wave in a density-jump medium to explain the unusual temporal 
feature of the GRB 030226 afterglow. A large density jump can produce a relativistic reverse 
shock, which leads naturally to (1) effective conversion of the initial kinetic energy to thermal 
energy and (2) rapid deceleration of the jet. Therefore, the assumption of one large density jump 
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has two advantages. Based on them, the observed rebrightening is understood as due to emissions 
from reverse and forward shocks forming in the interaction of the blast wave with the density 
jump, and the late-time afterglow is caused by sideways expansion of the jet. 

We thank the referee and X. Y. Wang for valuable comments. This work was supported by 
the National Natural Science Foundation of China (grant 10233010) and the National 973 Project 
(NKBRSF G19990754). 
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Fig. 1. — Fitting of the multiwavelength afterglow of GRB 030226 by assuming an ultrarelativistic 
blast wave which expands in a density-jump medium. The model parameters are seen in the 
text. Data are from Ando et al. 2003a, 2003b; Covino et al. 2003; Fatkhullin et al. 2003; 
Garnavich, von Braun & Stanek 2003; Greiner et al. 2003b ;Guarnieri et al. 2003; Maiorano et al. 
2003; Nysewander et al. 2003; Price & Warren 2003; Rumyantsev, Biryukov, & Pozanenko 2003; 
Rumyantsev, Sergeeva & Pozanenko 2003; Semkov 2003; von Braun, Garnavich & Stanek 2003a, 
2003b; Zeh et al. 2003. 



